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Abstract 

We discuss the ground state and the small-amplitude excitations of a degenerate vapour of 
fermionic atoms placed in two hyperfine states inside a spherical harmonic trap. An equations- 
of-motion approach is set up to discuss the hydrodynamic dissipation processes from the inter- 
actions between the two components of the fluid beyond mean-field theory and to emphasize 
analogies with spin dynamics and spin diffusion in a homogeneous Fermi liquid. The condi- 
tions for the establishment of a collisional regime via scattering against cold-atom impurities 
are analyzed. The equilibrium density profiles are then calculated for a two-component vapour 
of 40 K atoms: they are little modified by the interactions for presently relevant values of the 
system parameters, but spatial separation of the two components will spontaneously arise as 
the number of atoms in the trap is increased. The eigenmodes of collective oscillation in both 
the total particle number density and the concentration density are evaluated analytically in 
the special case of a symmetric two-component vapour in the collisional regime. The disper- 
sion relation of the surface modes for the total particle density reduces in this case to that of 
a one-component Fermi vapour, whereas the frequencies of all other modes are shifted by the 
interactions. 



PACS. 67.40.Db Quantum statistical theory; ground state, elementary excitations 



1 Introduction 

The experimental realization of Bose- Einstein condensation in confined vapours of alkali atoms 
jj], [], ^, has given impulse to the study of dilute quantal fluids, including vapours of fermionic 
atoms. Magneto-optical confinement of fermionic species has been reported for 6 Li || and 40 K 
|J. DeMarco et al. M have realized magnetic trapping of 40 K atoms in two different hyperfine 
states corresponding to \F = 9/2, = 9/2) and \F = 9/2, F z = 7/2), with the possibility 
of varying the relative ' concentration of these two components of the vapour up to selective 
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removal of one of them. Earlier experimental work on double Bose condensates [g, and some 
of the related theoretical work on the equilibrium state and on the excitation properties of 
bosonic mixtures jl^, jf^, [l3|, [TIJ [Ts| may also be recalled at this point. 

The s-wave collisions between pairs of fermions in the same hyperfine state are suppressed by 
the Pauli principle, so that to leading order only p-wave scattering and dipole-dipole magnetic 
interactions remain in a one-component, spin-polarized Fermi vapour. These effects are very 
weak at very low temperatures and the vapour may be treated as an ideal Fermi gas Jl6| , |l7| , [l8| 
|is| ]. In the two-component vapours studied by DeMarco et al. Q, however, s-wave scattering 
is operative between pairs of 40 K atoms in different hyperfine states. They have thus been 
able to measure the s-wave scattering length of 40 K, to observe directly the p-wave energy 
threshold law and to evaporatively cool the vapour down to 5 mK. While the s-wave scattering 
determined in this way for 40 K is repulsive (i.e. is described by a positive scattering length), 
a negative scattering length for 6 Li atoms holds promise of achieving a superfluid state in a 
mixture of 6 Li atoms prepared in two hyperfine states po| |. The insurgence of superfluidity 
may be revealed through the study of the elementary excitations of the vapour j^], 2^, . 

In the present work we extend to two-component interacting Fermi vapours in the normal 
(non-superfluid) state our former study of the small-amplitude excitationsof density fluctuations 
in an ideal Fermi gas confined in a harmonic trap at zero temperature [pifl. We make use of 
an equations-of-motion approach which is formulated in full generality in Sect. 2 in order to 
stress the analogies between the problem of present interest and that of spin dynamics and spin 
diffusion in a homogeneous Fermi liquid in a given state of partial spin polarization pt| . 
The nature of the assumptions which are adopted in our further calculations on confined Fermi 
vapours is made more explicit and justified by this discussion. 

In Sect. 3 we assume complete equilibrium for the two-component vapour under spherical 
harmonic confinement and evaluate the Thomas-Fermi ground-state densities upon relating the 
components of the kinetic stress tensor to the local densities by the ideal-gas formula. We give 
specific attention to three different cases, i.e. (i) the 40 K system studied by DeMarco et al. Yq, 
(ii) a strong-coupling regime in which the repulsive interactions between the two components 
of a symmetric vapour drive their spatial separation, and (iii) a simplified description of the 
weak-coupling regime in a symmetric vapour. By a symmetric vapour as treated in (ii) and (iii) 
we mean equal numbers of particles in the two components as well as equal masses and equal 
confinements, as is relevant in relation to the experiments of DeMarco et al. The form 

obtained in (iii) for the density profile is used in Sect. 4 to obtain an analytic determination of 
the eigenvectors and of the dispersion relation for both in-phase and out-of-phase oscillations 
in a symmetric vapour in the collisional regime. The role of the interactions in comparison 
with the vibrational properties of an ideal one-component Fermi gas is of main interest here. 
Finally, Sect. 5 gives a brief summary of our main results and offers some concluding remarks. 

2 Generalized quantum hydrodynamics in a two com- 
ponent Fermi fluid 

We review in this section some general properties of the dynamics of a two-component fluid with 
given equilibrium densities n a (r), a being a component index that we shall write as a = (|, J.) 
to stress the analogy with the problem of spin dynamics in a partially spin-polarized electron 
gas [Q. The Hamiltonian describing the fluid in the presence of external scalar potentials 
V„(r,t) is 
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where (r, t) are the field operators and </> CTlCT ' (r, r') the interatomic potentials. Redistributions 
of population in the two states are not allowed. 

The equations of motions for the partial particle densities n a (r, t) are obtained by a standard 
procedure (see e.g. pffl), involving (i) the derivation of the eq uati on of motion for the density 
matrix p CT (x, x';f) ==7\i£(r, i)i/> CT (r, t)} from the Hamiltonian ( |2.l[ ), and (ii) projection on the 
diagonal r = (x + x')/2. Setting r' = x — x', the result is 

m f n ^ t] = V( r, V< r) n^(r,t)+vW[n„(r,t)Vi r) Vf(r,f)] 

+ ^ydVvL r) {[vL r V^'(r,r')] (p.(r,t)p a ,(r',t)) c } , (2.2) 

a' 

where the convention of summation over repeated Cartesian indices in the derivatives has been 



adopted. In Eq. (2.2) we have defined the kinetic stress tensors 



n^(r, t) = - — vL r '' V< r,) p CT (r - r'/2, r + r'/2; t) I (2.3) 



and the mean-field potentials 



V?(r,t) = V !T (r,t) + ^ / d 3 r'<f>„ t „,(r,r')n a ,( r ',t) . (2.4) 



The non-mean-field effects are collected in the last term on the RHS of Eq. (2.2), where p a (r,t) 
is the density operator and (p CT (r, i)/9 CT /(r', t)) c is the cluster part of the density-density corre- 
lations. No assumption has as yet been made on the temperature of the fluid. 



The equilibrium equations determining the density profiles n CT (r) are obtained from Eq. (2.2) 
by taking the static limit. The equations of motion for the density fluctuations 5n CT (r, t) dr iven 
by weak external potentials are then obtained by writing n ff (r, t) = n (T (r) + 5n (T (r, t) in Eq. (2.2) 
and by linearizing it. We shall go through these steps in Section |§| and ^ for a dilute two- 
component Fermi vapour. Here we proceed to introduce the approximations that we shall 



make in the dynamical treatment of Sect, a by discussing the non-mean-field term in Eq. (2.2) 



2.1 Interdiffusion in the two component fluid 

We evaluate in this section the role of collisions between fluctuations in determining damp- 
ing of collective motions in the two-component fluid. As a preliminary we recall that the 
linearized equations of motion for the partial density fluctuations in the two-component fluid 
are conveniently transformed into those for the total particle density fluctuation 5n(v, t) and 
for the concentration fluctuation 5M(r,t) (the "magnetization" fluctuation in the electron gas 
analogue) by taking simple linear combinations of the two 8n& (r,t)'s (see e.g. fjjf ). 

We consider first the dynamics of small fluctuations around a homogeneous equilibrium 
state, where we can appeal to the treatment given by Caccamo et al. ]26| for the evaluation 
of the interdiffusion (or "spin diffusion") coefficient. Momentum conservation ensures that the 
only non-vanishing inverse relaxation time in the hydrodynamic limit is the interdiffusion one, 
say t^ m , which is written as 

t m\i = n-jM/irnn^ny) (2.5) 

where raj and rij. are the partial equilibrium densities, n — n-\ + % is the total den sity and we 
have assumed m = m\ = mj.. An exact expression for the quantity 7a/ in Eq. (J2.5T) is obtained 



from the non-mean-field term in Eq. (2.2) in the appropriate hydrodynamic limit. In the case 
of a central pair potential this reads 

y M = J d\ (k • V)fr t (r) a ^(R)^(R + r)) e (2 6) 

in a reference frame where the component f is at rest and the other component is flowing with 
a uniform drift velocity Wj.. 
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The Fourier transform of the non-equilibrium correlation function in Eq. ( |2.6|) is evaluated 
in a binary collision approximation by the decoupling procedure used by Baym [29] in treating 
the electrical resistance of metals (see also Kadanoff and Baym ]p0[). Namely, 



1 f°° A 

nW n (k) [S n (k ) a;)S u (-k J - W )-S TT (-k J - W )S u (k )W )] , (2.7) 
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where SWa^ka;) is the van Hove dynamic structure factor of each component in the non- 
equilibrium state. For the dilute Fermi fluid of present interest we can replace the interaction 



potential in Eq. (2.7) by a contact interaction and 5 , CT o-(ka;) by the ideal-gas value correspond- 



ing to a displaced Fermi sphere for the J. component. Following the lines of the calculation 



given in Ref. ]2(J and taking for simplicity n-\ — ny in Eq. (2.5), we find to leading order in the 
temperature T the result 

Tmm = (A7vma 2 n E 2 F /3h 3 )(T/T F ) 2 . (2.8) 

Here, an is the (triplet) scattering length, Ef is the Fermi energy and Tf = EfI^b- This 
result could also be obtained directly from Eq. (6.8) in Ref. pq] upon replacing a screened 
Coulomb interaction by a c onta ct interaction. 



The result given in Eq. (2.8) above for a homogeneous, two-component Fermi fluid can now 
be used for an estimate of the role of collisions in a confined Fermi fluid. We replace the Fermi 
energy Ef by its local value, which in the case of harmonic confinement in a spherical trap 
characterized by a frequency ujf is 

E F = (3N) 1/3 huj f (2.9) 

with N the total number of fermions. Hence, 

(w/tmm)" 1 = (4ir/3 1/3 )(N 1/3 a n /a ho ) 2 (T/T F ) 2 (2.10) 

a-ho = (h/mtuf) 1 / 2 being the harmonic-oscillator length. A similar result has been reported 
recently by Vichi and Stringari |n| from a collision-integral approach. 

In summary, because of momentum conservation the damping processes in the hydrody- 
namic limit of a two-component Fermi fluid are associated with collisions between the two 
components and affect only their relative motions. These processes vanish quadratically with 
decreasing temperature because of Fermi statistics (see also A collisional regime may 

nevertheless be established by scattering against impurities (see for instance the work of Ruck- 
enstein and Levy J3^] on spin dynamics in paramagnetic quantum fluids). We turn below to 
an estimate of these collision processes in the normal Fermi fluid of present interest. 



2.2 Collisional regime via impurity scattering 

A collisional regime is established in the low-temperature vapour for both in-phase and out-of- 
phase modes of motion of the two components when the inequality 

wr«l (2.11) 

holds, r being the collision time for scattering of fermions against impurities and ui being on 
the scale of the trap frequency (u ~ cof) for low- frequency modes. 

For an estimate of the needed number N s of scatterers we take the impurities as cold 
atoms with a mean velocity which is negligible relative to that of the fermions. We can then 
write t = l/v, I and v being the mean free path and the average speed of a fermion. We have 
I = E _1 , where E is the macroscopic cross-section given by E = n s a in terms of the density n s of 
scatterers and of the cross-section a for fermion-impurity scattering (see for instance Ref. 133]). 



4 



Setting n B = A r s (47ra| /3) 1 and a — 47ra 2 c with a sc the fermion-impurity scattering length, 
and taking v = (3Ef ■/4m) 1 ^ 2 with Ef given by Eq. (Eli), we find 



1/2 / N 2 



r _ (4/27) 
' iV s (3iV) 

For illustrative purposes we consider the case of 39 K or 41 K bosonic impurities in the gas of 
40 K fermions studied by DeMarco et al. Q (N ~ 10 7 and w/ ~ 209 s" 1 , the latter being the 
geometric mean of the radial and axial frequencies in the experiment). From the known values 
of the 39 K- 40 K and 41 K- 40 K scattering lengths (a B —3600 and a s ~ 93 Bohr radii, respectively) 
we find that a number of 39 K impurities of order N s ~ 10~ 6 iV, or of 41 K impurities of order 



N s ~ 10~ 3 iV, would suffice to verify the inequality (2.11) with u = uJf. 

We conclude, therefore, that a collisional regime can easily be established for the low- 
frequency excitations of trapped Fermi vapours. This regime reflects rather directly the quantal 
statistics of the vapour |m| and we study it for the two-component Fermi fluid in Sect. [l] below. 
Excitations in the collisionless regime are of less interest, since they mostly reflect the frequency 
of the trap |lj]. 



3 Equilibrium density profiles in spherical confine- 
ment 

As already discussed in Sect. [j], we treat a dilute two-component Fermi gas at zero tempera- 
ture in which only s-wave scattering between pairs of fermions in different hyperfine states is 
operative. This coupling is described by the parameter / = 4-Kh 2 afi/rn. In the experimentally 
relevant situation the two populations have not only the same mass but also essentially identical 
numbers and trap frequencies. However, we shall impose the equality ATf = JVj, = N/2 only 
later below. 

We take the gas as being statically in the equilibrium state and dynamically in the collisional 
regime. As in our earlier work |24j, we relate the kinetic stress tensor of each component to its 
local density by the homogeneous Fermi gas formula, 

n^(r,i) = S a0 ^A[n„(r,t)] 5/:i (3.1) 

where A = ft 2 (67r 2 ) 2/3 /2m. Such a local density approximation assumes that the length scale 
for the variation of the density profiles in space is large relative to the inverse Fermi wave 
number kj 1 and to the length c/w, with w the excitation frequency and c the velocity of sound 
propagation in the homogeneous fluid. 



The equilibrium density profiles are then easily obtained from the static limit of Eq. (2.2) 
in the mean-field approximation. They have the Thomas-Fermi form, 

n CT (r) = 6[e a - 14(r) - /Mr)] {A' 1 ^ - V a {v) - /«,(r)]} 3/2 , (3.2) 



where a denotes the component different from a. In Eq. (3.2) Vo-(r) are the static confining 
potentials and e a are the chemical potentials, to be determined from the condition N a = 
J d 3 rn CT (r). We emphasize that the JV^'s are fixed, i.e. these equations do not allow for 
redistributions of population in the two hyperfine states. 

In the same approximation the total energy of the vapour is the sum of three terms, i.e. a 
kinetic energy Ekin, a potential energy Eho and an interaction energy Ei„t- These are 



2 



E km = (6^ 2 ) 2 / 3 !^ / d 3 rK(r)] 5 / 3 , ( :!.:>,', 



E ho = ^2 I d 3 ™ CT (r)V CT (r) (3.4) 
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FIGURE 1: Density profiles n-\(r) and n^(r) (in units of aZ^) versus 
distance r from the centre of a spherical trap (in units of ato) in 
a mixture of 10 7 fermions at composition N^/N^ — 3. The long 
(short) dashed curves show the profiles corresponding to an s-wave 
scattering length o-fj. = 157 as (a^ = —157 as), relatively to the 
case of the ideal mixture (full curves). 



and 



Eim = J^J d 3 r-n T (r)n;(r) 



(3.5) 

Thes e wi ll be helpful in understanding the behaviour of the vapour at strong coupling in 



Sect. 3.2 



3.1 An illustrative example for a weakly coupled 40 K vapour 

Figure |l| reports the numerical results that we obtain from Eq. ([^) for the density profiles 
in a gas subject to spherical harmonic confinement, with system parameters chosen after the 
experiment of DeMarco et al. @ (w/ = (w||wi) 1/3 = 209 s" 1 , N = 10 7 and a n = 157 Bohr 
radii) but at composition Nf/Ni = 3. The cases a-fj. = —157 as and Oft = are also shown. 

Evidently, the effects of the interactions are small in this situation and very simple to 
understand: a repulsion (attraction) between the two components disfavours (favours) their 
overlap in the central part of the trap. 

3.2 Spatial separation of the two components at strong coupling 

We may expect that in the case of repulsive interactions, with increasing coupling stre ngt h 



and still barring transitions between the two hyperfine levels as already noted under Eq. (3.2), 
the gas will be led to diminish its total energy by reducing the spatial overlap between the 
two components. Figure |^ shows how such symmetry breaking occurs in a spherical trap, for 
JVf = iVj and under the condition that overall spherical symmetry be maintained. In this case 
one component is pushed away from the centre of the trap and the gas configuration becomes 
that of a central core enriched in one component and surrounded by a spherical shell enriched 
in the other. 
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Figure 2: Density profiles (in units of a h ^) versus distance r from 
the centre of a spherical trap (in units of a/j D ) in a symmetric mixture 
of fermions, at various values of the coupling strength parameter 
7. In (a), the profiles of the two components are still in complete 
overlap. Spatial symmetry breaking is first visible in (b). Separation 
of the two components (shown by full and dashed curves) continues 
through (c) and (d). 
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Figure 3: The ratio of the interaction energy Ei nt to the sum 
Eho + Ekin of the harmonic-oscillator and kinetic energies, plotted 
against the coupling strength 7 ~ 0.5A^ 1 / 6 (a||/a/ lo ) for a symmetric 
mixture of fermions in spherical confinement. The inset shows the 
first derivative of the same function. 



The symmetry breaking is driven by the competition between the repulsive interaction 
energy, favouring spatial separation of the components, and the kinetic energy disfavouring 
localization. Evidently and in contrast to the behaviour illustrated in Figure |2|, a gas confined 
in an axially symmetric trap will tend to reduce its energy via relative shifts in the centres 
of the two clouds. A rich phase diagram will ensue if "spin flips" between the two hyperfine 
states are also allowed. In the following we estimate the critical coupling strength at which 
spatial symmetry breaking occurs in terms of the number of fermions (or alternatively of the 
Of 4, scattering length) for the case illustrated in Figure bl 

The four cases of density profiles illustrated in Figure Tu are labelled by a parameter 7, which 
is defined by 

1 = aN 1/(i (a n /a ho ) (3.6) 

with a = 2 1/3 3 1/6 (8192/2835tt 2 ) Q. In fact the value of 7 in Eq. is obtained as 7 = 

[Eint I (Eho + Ekin)]o when the ratio Ei nt / (Eho + Ekin) is calculated from the density profile of 
the Fermi gas in the absence of interactions. Figure Hreports the true values of Eint / (Eho + 
Ekin) against 7 for the system described in FigureEi The (obvious) linear shape of this 
function at weak couplings gently bends over with increasing coupling, until an almost sharp 
break occurs at spatial symmetry breaking. This is emphasized in the inset in Figure ^ giving 
the derivative of Eint I (Eho + Ekin) with respect to 7. We have checked that the same plot is 
obtained for Eint / '(Eho + Ekin) by varying a^i/aho at constant N and by varying N at constant 

a 7i/ a ho- 

In the experiments of DeMarco et al. [jij on 40 K with N ~ 10 7 , the value of 7 is 7 ~ 0.022 
i.e. still very far from the critical value 7 C ~ 0.535 for the symmetry breaking illustrated in 
Figures |^ and ^. The weak dependence of 7 on N in Eq. (3.6) implies that a number of 40 K 
atoms of order 10 15 would have to be reached if all other system parameters remain the same. 
A parallel increase in the ratio an/aho as suggested by Eq. (3.6) would evidently be helpful in 
relaxing such stringent condition. 
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Figure 4: Illustrating the accuracy of the approximate form ( |3.8| ) 
of the particle distribution \-Kr 2 n(r) (dashed curves) relative to the 
full Thomas- Fermi profile (full curves), for a symmetric mixture of 
fcrmions in the two cases 7 = ±0.022. 



3.3 Approximate form of the density profile for a symmetric 
vapour at weak coupling 

In the case iVj = iVj. the shape of the total density profile n(r) = nj (r) +nj. (r) at weak coupling 
is well represented by a form which is suitable for t he analytic study of the eigenmodes of the 
gas that we report in Sect. ^ below. From Eq. (3.2) the Thomas- Fermi density profile is 



i(r) = 2A 



-3/2 



1 2 2 

e F - -mu f r 



- -n(r) 
2 v ; 



3/2 



(3.7) 



for r < Rf, where Rf = (2e F /rnuj 2 f ) 1 |/2 a nd e F = £f = ej is the chemical potential of the 



mixture. At weak coupling the profile (3.7) can be approximated by 



i(r) = 



_87V_ 

7V 2 R F 



(1- 



r 2 /R 2 F ) 3/2 t 



(l-r 2 /R%) 



(3.8) 



where and st ill a re the true chemical potential and the Fermi radius in the interacting 
mixture. The form (3.8) is adjusted to preserv e no rmalization to N as well as the value of Rf - 
Figure [| compares the approximate form (3.8) with the correct Thomas- Fermi form (3.7) 
at 7 = ±0.022 by plotting the function 4irr 2 n(r). The phase-space factor 4nr 2 masks the small 
differences that would be present in the two forms of n(r) near the centre of the trap. On the 
other hand, preserving the correct value of the Fermi radius in the approximate form (3.8) is 



crucial in view of the boundary conditions to be imposed in the determination of the eigenmodes 
of the vapour. 



As a final remark we notice that the profile in Eq. (3J3) has the same form as for an ideal 
one-component Fermi gas j34|. This fact is crucial for the analytic treatment of the dynamics 
of density fluctuations in a weakly coupled symmetric mixture, that we give in the next section. 
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4 Dynamics of density fluctuations 



The equations of motions ( p.2[ ), after linearization in the partial density fluctuations and adopt- 
ing (i) the mean-field approximation ((p tT ( r, t) p a i (r',t)) c ~ 0) and (ii) the local density approx- 
imation for the kinetic stress tensor (Eq. (3.f), reduce to 



3V(r) 



^An 2 J 3 (v)5n a (v,t) 



= V" 

+ V ■ {<5n CT (r, t) V[V CT (r) + fn g (r, t)] + fn a (r)V6n g (r, t)} 



(4.1) 



at resonance (i.e. for 14(r,£) — V CT (r)). With the help of the equilibrium conditions (3.2), and 
taking Fourier transforms with respect to the time variable, Eq. (^^) can be written as 



mui 2 Sria- (r , lu) 
+ {[e S 



\A [2n 2 J 3 (v)V 2 + V{n 2 J 3 ) ■ V - V 2 {n 2 J 3 )] 5n a (v,uj) 



3 

14 (r 



An 2/3 W 2 



■V[V g (r)+An 2/3 ]-V 



7 }Sn & (r,Lu) (4.2) 



Evidently, Eq. (4.2) describes a two-by-two eigenvalue problem for the coupled partial density 
fluctuations, which is to be solved numerically in the general case. 

The problem is considerably simplified in the case of a symmetric mixture (rri| = mj, 
ATf = N± and V\ (r) = Vy (r), these conditions be ing well satisfied in the experiments of DeMarco 
et al. Q) . In this case the dynamical equations ( |4.2| ) lead to separate equations of motion for the 
total density fluctuations Sn(r, lu) = <5rif (r, ui) +<5rij(r, u>) and for the concentration fluctuations 
8M(v,oj) — Sni(r,uj) — 8n±(v,oj). The eigenvalue equation for 8n(v,u)) reads 



mui 5n(r,ui) 



7 ■ (^A(n/2) 2/3 VSn(r,to) ~ ^{V(n/2) 2/3 )8r, 
+ (e F - V(r) - A(n/2) 2/3 )VSn(r, u))) 
Similarly, the eigenvalue equation for SM(r, ui) is 



(4.3) 



mui 5M(r, uj) 



V- 



^A{n/2) 2/3 X75M 



A 



(V(n/2) 2/3 )8M(r,iu) 



(e F - V(r) - A{n/2) 2 ' 3 ) VSM (r, w)) 



(4.4) 



Equations (f.2) and (4.4) can be solved analytically by the technique used in our earlier work on 
the ideal one-component Fermi gas pi| , if the form (i.8) is adopted for the equilibrium profile. 
As already discussed in Sect. t^Jj, Eq. (3.8) becomes accurate at small coupling. We shall 



again impose that the solutions vanish continuously at the cloud boundary, as a consequence of 
Fermi statistics giving a high cost in kinetic energy to rapid variations of the densities in space. 
For both in-phase and out-of-phase motions of the two components, the frequency eigenvalues 
depend on a parameter C given by 



C = (3N) 2/3 (hco f /€ F ) 



(4.5) 



This quantity is the square of the ratio of the ideal Fermi energy to the true Fermi energy and 
hence, in the case of repulsive interactions where e F increases with the scattering length, is 
limite d fr om a bove by the inequality C < 1. All the mathematical details of the solution of 
Eqs. (4.3) and (4.4) are given in Appendix Here we report only the main results. 



4.1 Small oscillations of total density fluctuations 

The eigenfunctions of the total density oscillations vanish at the Fermi radius r = R F provided 
C < 3, in a wa y wh ich depends on the parameter C and hence on the strength of the interactions 
(see Appendix A.l for their detailed expressions). 
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The corresponding eigenfrequencies are labelled by the angular momentum number I and 
by an integer n representing the number of internal nodes in the density fluctuation profile. 
The dispersion relation is 

(wni/w/) 2 = l + 2n+ ~(3 - C)(2n + 21 + 1) (4.6) 



The ideal Fermi gas limit corresponds to C = f and in t his case Eq. (4.6) yields back our earlier 
result [p4| . More generally, the dispersion relation (4.6) reduces to that of the ideal Fermi gas 



only for the surface modes (i.e. for n = 0). Instead, the frequencies of the modes with n > 
are shifted by the interactions. 

The above results are easily extended to evaluate the low-frequency modes in an axially 
symmetric confinement (see e.g. As discussed in Sect. ^, damping of these modes will set 

in when ujt ~ 1 for scattering against cold-atom impurities. 

4.2 Small oscillations of concentration fluctuations 



The eigenfunctions vanish at the Fermi radius only if C > 3/5 (see Appendix A. 2 for the details). 
This inequality marks the breakdown of the present approximation in the case of repulsive 
interactions. Under this restriction the dispersion relation for concentration fluctuations having 
n internal nodes is 

{u n i/uf) 2 = (l + 2n)(2C - 1) + -(5(7 - 3)(2n + 21 + 1) (4.7) 

Of course, for both surface (n = 0) and bulk (n ^ 0) modes these frequencies differ from those 
of density fluctuations in the ideal one-component Fermi gas. 

Damping of these modes will arise not o nly from scattering against cold-atom impurities 



but also from thermal excitations (see Sect. 2.1. A detailed discussion of the latter damping 



for the spin dipole excitation has been given by Vichi and Stringari 



5 Summary and concluding remarks 

In summary, the focus of this work has been on two-component mixtures of fermionic atoms in 
dilute-vapour states and subject to spherical harmonic confinement at zero temperature. The 
main motivation has come from the experiments of DeMarco et al. [6] on vapours of 40 K atoms 
magnetically trapped in two different hyperfine states. 

The generalized hydrodynamic equations of the mixture have allowed us to discuss the 
damping mechanisms from correlations between partial density fluctuations beyond mean field 
terms. Because of momentum conservation in a pure two-component Fermi fluid the dissipation 
processes in the hydrodynamic limit are associated with collisions between the two components 
and these vanish quadratically with temperature because of Fermi statistics. We have then 
discussed how a collisional regime may nevertheless arise for both global and relative density 
fluctuations at very low temperature from collisions of the Fermi fluid against cold impurity 
atoms. We have seen that the establishment of a collisional regime, in which the dynamical be- 
haviour of the fluid reflects the quantal statistics, is not subject to especially severe restrictions 
on the strength of the fermion-impurity scattering nor on the number of impurities. 

We have then evaluated ground-state properties and small-amplitude excitations of such a 
two-component Fermi fluid in a collisional regime. We have shown that, whereas the role of 
the interactions in determining the equilibrium density profiles is still very weak in the cases 
experimentally studied so far, a rich phase diagram will emerge as the coupling strength is 
increased and/or redistribution of the components between magnetic states becomes allowed. 
The relevant coupling strength depends in a simple manner on the number of fermions in the 
trap and on the ratio of the scattering length to the harmonic-oscillator length. Finally, we have 
shown how the problem of small-amplitude oscillations of both the total particle density and 
the concentration density in a weakly coupled symmetric mixture is amenable to full analytic 
solution in parallel with the analogous problem for an ideal Fermi gas. 
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A Solution of equations (||73|) and (gTg ) 



We give in this Appendix the details of the analytic solution of the eigenvalue equations (4.3) 
and (1.4) and the expressions for their eigenfunctions. 



A.l Total density fluctuations 

From Eqs. (|J) and (|J) we find 

[&{L0/i0ff +CV^(i 2 )]fe(x,[j) 



+(3 - C)(l - x 2 )X7 2 Jn(x, oj) - (3 - 2C)V(x ) • V<5n(x, w) = 



(A.l) 



where x = r/R F . The solutions of Eq. (A.l) have the form <5n(x,w) = x l F(x 2 )Y, m ( 8 <j)) 
because of spherical symmetry. Setting x 2 — y, we determine the function F(y) from Eq. (A.l) 
by means of the Fuchs method for solving an ordinary differential equation in a series form 
around regular singular points pah. This method sets 



F(») = (1-V)*X>*(1-V)* 



(A.2) 



and yields s — C/(3 — C) together with the recurrence relation for the coefficients at, 



2(s + fc + l)(s + fc-36 + 2) 



a-k+i 



-3(36 - 1) + 1(2 - 36) + (s + k)[2{s + k - 1) + 21 - 66 + 7] 



(A.3) 



Here, 6 = (3 — C) 1 . The eigenfunctions vanish at th e boundary for C < 3. 

The eigenfrequencies are obtained from Eq. (A.3) by asking that the solutions reduce to 



polynomials of degree n + s, i.e. a n +i = for an integer n representing the number of inte rna l 
nodes of the density fluctuation profile. This yields the dispersion relation reported in Eq. (4.6) 
of the main text. 



A.2 Concentration fluctuations 



From Eqs. (4.4) and (3.8) we get 

[6(w/w/) 2 + C\7 2 x (x 2 )}5M(x, to) 

+(5C - 3)(1 - x 2 )V 2 JM(x,uj) + (3 - 4C)V(a; 2 ) • V<5M(x,w) = (A.4) 

We look for solutions having the form SM(x,lu) — x l G(x 2 )Y l m (8, cj>) and set x 2 = y to find the 
differential equation obeyed by the function G(y), 

2(5C - %(1 - + [(5C - 3)(2J + 3)(1 - y) + 2y(3 ~ AC)] dG ' 



dy 2 



dy 



+ (3{ij/LUf) 2 + 3C + 1(3 - AC)) G(y) = (A.5) 
Following again the Fuchs method we set 

G( y ) = (i- y yY d a ^ 1 -y) k (A.6) 



and from the indicial equation for Eq. ( A5) we find s = C/(5C — 3). Therefore, the solutions 
will vanish at the boundary of the cloud only if C > 3/5. 
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The coefficients of the series in Eq. (A.6) obey the recurrence reiation 



2(s + fc + l)[(s + fc)(5C-3) -3 + 4C]^±i = -3(uj/LUf) 2 

ak 

-[3C + 1(3 - AC)] + (s + k)[(5C - 3) (2s + 2k + 21 + 1) + 2(4C - 3)] . (A.7) 

By asking again for polynomial solutions, we obtain the dispersion relation for concentration 
fluctuations having n internal nodes as given in Eq. J4.7| ) of the main text. 

References 

M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Science 
269, 198 (1995). 

K. B. Davis, M. O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Durfee, D. M. Kurn, 
and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995). 

C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. Hulet, Phys. Rev. Lett. 75, 1687 
(1995), see also C. C. Bradley, C. A. Sackett, and R. G. Hulet, Phys. Rev. Lett. 78, 985 
(1997). 

C. C. Bradley, C. A. Sackett, and R. G. Hulet, Phys. Rev. Lett. 78, 985 (1997). 

W. Mc Alexander, E. Abraham, N. Ritchie, C. J. Williams, H. Stoof, and R. Hulet, Phys. 
Rev. A 51, R871 (1995). 

F. S. Cataliotti, E. A. Cornell, C. Fort, M. Inguscio, F. Marin, M. Prevedelli, L. Ricci, 
and G. M. Tino, Phys. Rev. A 57, 1136 (1998). 

B. DeMarco, J. L. Bohm, J. P. Burke, Jr., M. Holland, and D. S. Jin, Phys. Rev. Lett. 
82, 4208 (1999). 

C. J. Myatt, E. A. Burt, R. W. Christ, E. A. Cornell, and C. E. Wieman, Phys. Rev. Lett 
78, 586 (1997). 

D. S. Hall, M. R. Matthews, J. R. Ensher, C. E. Wieman, and E. A. Cornell, Phys. Rev. 
Lett. 81, 1539 (1998). 

B. Esry, C. Greene, J. B. Jr., and J. Bohn, Phys. Rev. Lett. 78, 3594 (1997). 
R. Graham and D. Walls, Phys. Rev. A 57, 484 (1998). 
B. Esry and C. Greene, Phys. Rev. A 57, 1265 (1998). 
P. Ohberg and S. Stenholm, Phys. Rev. A 57, 1272 (1998). 
H. Pu and N. P. Bigelow, Phys. Rev. Lett. 80, 1130 (1998). 
H. Pu and N. Bigelow, Phys. Rev. Lett. 80, 1134 (1998). 
D. A. Butts and D. S. Rokhsar, Phys. Rev. A 55, 4346 (1997). 
J. Schneider and H. Wallis, Phys. Rev. A 57, 1253 (1998). 



H. T. C. Stoof and M. Houbiers, cond-mat/9808171 (unpublished) 



M. Amoruso, A. Minguzzi, S. Stringari, M. P. Tosi, and L.Vichi, Eur. Phys. J. D 4, 261 
(1998). 

M. Houbiers, R. Ferweda, H. T. C. Stoof, W. I. McAlexander, C. A. Sackett, and R. G. 
Hulet, Phys. Rev. A 56, 4864 (1997). 

W. Zhang, C. A. Sackett, and R. G. Hulet, Phys. Rev. A 60, 504 (1999). 



M. A. Baranov and D. S. Petrov, ;ond-mat/9901108 (unpublished). 



G. Bruun and C. W. Clark, cond-mat/9906392 (unpublished) 



M. Amoruso, I. Meccoli, A. Minguzzi, and M. Tosi, |cond-mat/9907370| and Eur. Phys. J. 
D, in press (unpublished). 



13 



A. Abrikosov and I. Khalatnikov, Rep. Progr. Phys. 22, 329 (1959). 

C. Caccamo, G. Pizzimenti, and M. Tosi, N. Cimento B 31, 329 (1976). 
A. Minguzzi, M. Chiofalo, and M. Tosi, Phys. Lett. A 236, 237 (1997). 

D. Kim, H. Praddaude, and B. Schwartz, Phys. Rev. Lett. 23, 419 (1969). 
G. Baym, Phys. Rev. 135, 1691 (1964). 

L. Kadanoff and G. Baym, Quantum Statistical Mechanics (Benjamin, New York, 1962). 
L. Vichi and S. Stringari, fcond-mat/9905154 and Phys. Rev. A, in press (unpublished). 



A. Ruckenstein and L. Levy, Phys. Rev. B 39, 183 (1989). 

A. Weinberg and E. Wigner, The Physical Theory of Neutron Chain Reactors (University 
of Chicago Press, Chicago, 1958), p. 24. 

K. M0lmer, Phys. Rev. Lett. 80, 1804 (1998). 

C. M. Bender and S. A. Orszag, Advanced Mathematical Methods for Scientists and Engi- 
neers (McGraw-Hill, New York, 1978). 



14 



